Four-nucleon scattering: ab initio calculations in momentum space 
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The four-body equations of Alt, Grassberger and Sandhas are solved for n-'^H scattering at en- 
ergies below three-body breakup threshold using various realistic interactions including one derived 
from chiral perturbation theory. After partial wave decomposition the equations are three-variable 
integral equations that are solved numerically without any approximations beyond the usual dis- 
cretization of continuum variables on a finite momentum mesh. Large number of two-, three- and 
four-nucleon partial waves are considered until the convergence of the results is obtained. The total 
n-"^H cross section data in the resonance region is not described by the calculations which confirms 
previous findings by other groups. Nevertheless the numbers we get are slightly higher and closer to 
the data than previously found and depend on the choice of the two-nucleon potential. Correlations 
between the Ay deficiency in n-d elastic scattering and the total n-^H cross section are studied. 

PACS numbers: 21.45.-(-v, 21.30.-x, 24.70.-(-s, 25.10.+S 



I. INTRODUCTION 

The four-nucleon {AN) scattering problem gives rise 
to the simplest set of nuclear reactions that shows the 
complexity of heavier systems. The neutron-'^H [n-^H) 
and proton-'^He (p-^He) scattering is dominated by the 
total isospin T — 1 states while deuteron-deuteron (d-d) 
scattering by the T — Q states; the reactions n-'^He and 
p-^H involve both T = and T = 1 and are coupled 
to d-d in T = 0. Due to the charge dependence of the 
hadronic and electromagnetic interaction a small admix- 
ture of T = 2 states is also present. In AN scattering 
the Coulomb interaction is paramount not only to treat 
p-'^He but also to separate the n-'^He threshold from p-^H 
and at the same time avoid a second excited state of the 
a particle a few keV bellow the lowest scattering thresh- 
old. All these complex features make the AN scattering 
problem not only a natural theoretical laboratory to test 
different force models of the nuclear interaction, but also 
the next step in the pursuit of very accurate ah initio cal- 
culations of the iV-body scattering problem after the ex- 
tensive work on the three-nucleon (3iV) system that has 
takenplace in the past twenty years by several groups 

In Refs. [3, Q all the reactions mentioned above were 
studied in the framework of Alt, Grassberger, and Sand- 
has (AGS) equations Q using the rank one representa- 
tion of realistic two-nucleon (2A^) force models together 
with a high rank representation of all 3iV subsystem am- 
plitudes; the Coulomb interaction was neglected. This 
led to one-variable integral equations whose predictive 
power was limited to the quality of the involved ap- 
proximations. The calculations showed large discrep- 
ancies with data, namely nucleon analyzing power Ay 

in n-^He scattering, tensor observables in ^II((i, n)'^IIe 
and ^H((i, d)^Yi and the differential cross section for 
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^II((i, n)'^IIe, but one surprising success in describing the 
total cross section at for n-'^H scattering in the reso- 
nance region where at neutron lab energy En — 3.5 
MeV at rises to about 2.45 b Calculations by the 
Grenoble group Q using coordinate-space solutions of 
the Faddeev-YakulDOvsky equations 0] showed, on the 
contrary, that realistic interactions missed the total cross 
section peak by at least 0.2 b. Although these calcula- 
tions carried out no approximation on the treatment of 
the 2A^ interaction, they were limited vis-a-vis Ref. Q on 
the number of ?>N and 4A'' partial waves. 

Although the issue was recently clarified [l^, [ll| by 
comparing to an independent calculation by the Pisa 
group that uses the Kohn variational method, together 
with hyperspherical harmonics, further studies based on 
the AGS equations are needed to settle this important 
problem because some of the results by the Grenoble and 
Pisa groups may be still of limited accuracy given the 
number of included 27V, 2>N , and AN partial waves. Fur- 
ther investigations are also needed for the understanding 
of other 4A^ reactions such as p-'^He, n-'^He and d-d where 
large discrepancies with data were previously found. One 
fundamental issue underlying four-nucleon physics is the 
existence of correlations between 3A^ and 4A^ observables. 
One of the best known is the Tjon line [l^ which cor- 
relates the binding energies of '^H with *He; another one 
involves the triton binding energy and the singlet (triplet) 
n-^M scattering length [13!|. Nevertheless, other correla- 
tions may exist: one could ask if the persistent Ay prob- 
lem in n-d scattering is in any way related to the failure 
to reproduce at in n-'^H scattering in the resonance re- 
gion, or to the Ay problem in p-^He [14]; does resolving 
the former also solves the latter? 

Therefore we present here a new numerical approach 
to the solution of the AGS equations that is both nu- 
merically exact and extremely fast in terms of CPU-time 
demand. Since the 2N transition matrix (t-matrix) is 
treated exactly, the equations we solve are, after partial 
wave decomposition, three- variable integral equations. 
The three Jacobi momentum variables in 1 + 3 and 2-1-2 
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configurations are discretized on a finite mesh and the 
number of 2N , 3N and AN partial waves increased up 
to what is needed for the full convergence of the observ- 
ables. The present approach also allows for the inclusion 
of charge-dependent interactions as well as A degrees of 
freedom that lead to an effective 3A^ force. Furthermore, 
using the method recently proposed to treat the Coulomb 
force in p-d elastic scattering and breakup [l^ [l^, [l^ , 
we have already obtained preliminary results for p-^He 
elastic scattering observables jl8j] with the Coulomb po- 
tential between the three protons included. 

In Sec.|TT]we discuss the integral equations we solve for 
n-^H scattering and in Sec. lIIII we show the results of our 
most complete calculations, leaving tests of benchmark 
to the Appendix. Finally conclusions come in Sec. IIVI 

II. EQUATIONS 

As initially proposed by Alt, Grassberger and Sandhas 
[sl and later reviewed for the purpose of practical appli- 
cations in Ref. 19], the four-particle scattering equations 
may be written in a matrix form 

U = V + VgoU, (la) 
W|<I>p„)=V|*p„>, (lb) 

l*Po> = l*Po>+^oV|*p„), (IC) 

where \^pg) is the initial channel state, jvl'po) the full 
scattering state, and po defines the two-body entrance 
channel. Both of them have 18 components, and the 
transition operator U as well as V and Qq are 18 x 18 
matrix operators with components 

[Kf = iGoUGo)-^S^pS,,, (2a) 
[^o],^'' = Go t. Go U^^Go t, Go 6„p. (2b) 

As usual, a{p) denotes two-cluster partitions of 1 -f 3 
or 2 -t- 2 type and the pair interactions. Go is the 
four free particle Green's function, ti is the two-particle 
t-matrix embedded in four-particle space. Sap = 1 — 5a p, 
and U^j are the subsystem transition operators 

C/- = 1 % + ^'^ Go U^^ (3) 

fe 

of 1 -t- 3 or 2 -f 2 type, depending on ct. If ct is a 1 + 3 parti- 
tion, Ufj corresponds to the usual AGS transition matrix 
for the three interacting particles that are internal to a. 
For a of 2-1-2 type lH^j does not correspond to any physical 
process. The components of the initial/final two-cluster 
states [|$po)]f = |0f°)^ppo are the Faddeev components 
of the cluster bound state wave function times a plane 
wave of momentum Pp^ between clusters whose depen- 
dence is suppressed in our notation, 

m^GouY.~^r.^m■ (4) 

k 



The great advantage of AGS equations over the 
Yakubovsky equations is that on-shell matrix elements of 
U between initial |$po) and final \'^ao) states with rela- 
tive two-cluster momenta Pp^ and pj^,, lead automatically 
to the corresponding scattering amplitudes 

(p;jT-'"'«|PpJ = ($.jZ^|$pJ (5a) 

= Y.^^TPT"\^'°)- (5b) 

For four identical particles the AGS equations reduce 
to 2 X 2 matrix equations since there are only two distinct 
partitions, one of 1 -I- 3 type and one of 2 -f 2 type, which 
we choose to be (12,3)4 and (12) (34); in the following we 
denote them by a = 1 and a = 2, respectively. In this 
case the equations may be conveniently written using the 
permutation operators Pab of particles a and b as it was 
done first in Refs. [13, [2l[ for the four-nucleon bound 
state. After the symmetrization of the four-nuclcon scat- 
tering equations ([T]) we obtain equations of the same form 
but with new definitions for the symmetrized operators 



V"^-(GotGo)-i(l„/3- 5^1^34), (6a) 

^GotGoU'^GotGoS^p. (6b) 

Here t is the pair (12) t-matrix, [/" the symmetrized 1+3 
or 2-|-2 subsystem transition operators 

t/" -P„Goi+F„tGoC/", (7) 

and Pa the permutation operators given by 

Pi = P = Pi2 P23 + Pi3 P23, (8a) 

P2=P = Pl3P2i- (8b) 



The basis states are antisymmetric under exchange of 
two particles in subsystem (12) for 1 + 3 partition and 
in (12) and (34) for 2 + 2 partition. The symmetrized 
initial/final two-cluster state components are 

|0/') = GotP;3|0O- (9) 

The scattering amplitudes are obtained as 

{p'jT'^^\pp) = Sap{r\U''^\cl^^), (10) 

where Sap is a symmetrization factor; Saa = Sa is the 
number of pairs internal to the partition a, i.e., 6*1 = 3 
and 5*2 = 2, and 5i2 = 2521 = 2a/3. 

Since the present paper is confined to n-^H scattering, 
we write down explicitly only the equations for the 1 + 
3^1 + 3 and 1 + 3^2 + 2 transition operators 

- {GotGo)-'P34- P34U'GotGoU'^ 

+ U^GotGoU^\ (11a) 

U^^ = {GotGo)-\l - P3i) + {I - P3i)U^GotGoU^\ 

(lib) 

The equations coupling U^"^ and U"^^ share an identical 
kernel but different inhomogeneous terms. 
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FIG. 1: (Color online) Angular momentum quantum num- 
bers for 1 + 3 and 2 + 2 basis states. 



After the partial wave expansion Eqs. pT|l form a 
set of coupled integral equations with three variables 
corresponding to the Jacobi momenta k^, ky and k^; 
the associated orbital angular momenta are denoted 
by Ixi E^nd Iz, respectively. They are depicted in 
Fig. [T] for 1 + 3 and 2 + 2 configurations together 
with the pair total angular momentum / and the 
the three-particle subsystem total angular momentum 
J. The states of total angular momentum J7 are de- 
fined as \kxkykz[h{{ly[{lxSx)I Sy]Sy}Jsz)Sz]JM) 
for the 1 + 3 configuration and 
\kxkykz{lz{{lxSx)I[ly{sySz)Sy]r}Sz)JM) foT the 
2 + 2, where Sy and Sz are the spins of nucleons 3 and 
4, and Sx, Sy, and Sz are channel spins of two-, three-, 
and four-particle system. In all calculations / and /' run 
over the same set of quantum numbers. 

By the discretization of the momentum variables the 
integral equations may be turned into a system of linear 
equations but the direct solution is not possible because 
of the huge dimension. Therefore, in close analogy with 
three-nucleon scattering, we calculate the Neumann se- 
ries for the on-shell matrix elements of the transition op- 
erators dTTI) and sum by the Pade method [l^] . The Pade 
summation algorithm we use is described in Ref . ^23| . We 
work with the half-shell transition operators in the form 



(12) 



such that the on-shell elements are {(j)'^\l4"^\(j)^) = 
{^a\Xai3) with the auxihary states \^a) = Go^l^") ~ 
tP„|</>")- Defining Q„ = GnW^Got and using Eq. ^ 
for the inhomogeneous terms in order to eliminate 
(GotG'o)~\ Eqs. HI]) become 

l^ll) = - PsiPlGol^l) - P34gi|^ll) + Q2\X2l), 

(13a) 

1^21) = (1 - P34)PlGo|a) + (1 - -P34)Ql|^ll). (13b) 



In practical calculations, in order to accelerate the con- 
vergence of the Pade summation, it is advantageous to 
substitute Eq. (|13bp into Eq. p3ap yielding the Neumann 



l^i?^> = (l-^34)PlGo|6), 



n=0 



(14a) 

(14b) 
(14c) 

xi^^) = il-Ps,)Qi\xir'), (14d) 

l^li') = - P34Qi\xir'^) + Q2\X^'i^), (14e) 

which requires 1 + 3 and 2 + 2 subsystem transition op- 
erators [/", contained in Qa, fully off-shell at different 
energies. Explicit calculation of U°' is not only very time 
consuming but also requires large storage devices. There- 
fore, except at the bound state poles, we do not calcu- 
late the full off-shell transition matrices explicitly. 
Instead, we rewrite Eq. ([7]) as a Neumann series 



(15) 



resulting a corresponding Neumann series for the solution 
vectors in Eqs. (HH), i.e.. 



1^, 

\x 



(n,0)> 
Q/3 / 

(",'")\ 
a/3 / 



\X 



1 

a/3 / 



PaGo t \X 



(ri,r-l)> 
a/3 / 



(16a) 

(16b) 
(16c) 



where the summation again has to be performed using 
the Pade method. Usually, 6 to 18 Pade iteration steps 
are required for the convergence in Eqs. - pB]) . At 
the bound state poles the subsystem transition operators 



U"=Pc.\U, 



Set 



E + iO 



E, 



-{^a\Pa, (17) 



where E is the available four-nucleon energy, Eb the 
binding energy, and Hq the kinetic energy operator for 
the relative motion of the two clusters. 

Thus, compared to the calculation of full off-shell [/", 
the method we are using avoids storage problems and 
also significantly reduces the number of required floating 
point operations, since it is essentially a calculation of 
half-shell matrix elements for a number of driving terms 
that are considerably fewer than the linear dimension 
of the discretized U". A further advantage is that the 
matrices corresponding to the operators Pa, Gq and t in 
Eq. (|16cp have block-diagonal structure whereas U" is a 
full matrix. 

The calculation of the Neumann series ([T6|) for a = 1 
is what we are doing in three-nucleon scattering and is 
described in great detail in Refs. [H,!!^. The specific 
representation of the permutation operator Pi where the 
initial and final state momenta ky are chosen as indepen- 
dent variables requires the interpolation in the momen- 
tum kx for the quantities on both sides of PiGq, i.e., for 



4 



t or (^Q,|. Two interpolation methods using Chebyshev 
polynomials and spline functions were used in Ref. [2^ : 
in the context of four-nucleon equations where one has 
to work with 1 + 3 and 2 + 2 basis states the spline in- 
terpolation is more convenient. 

The calculation of the Neumann series p6|) for a = 2 
is straightforward because of the very simple form of the 
permutation operator P2. 

Finally, the application of the permutation operator 



P-iA as well as the transformation of \X 



(»)\ 



from 1 



34 aa well ao uiic ,.iciiioiL.iiiiauiuii kjl k^^^ , 

basis to 2 + 2 or vice versa has a structure similar to 
that of Pi, resulting in a similar treatment. The specific 
representation of P34, i.e., 

{kxkykz\P3,\Kk'yk'z) = ^-^^^^^^ dy P3,{kz,k'z,y) 

^{^v ~ kyjkz, kz,y)) 
_ ^ 

j(fc;-fc;(fc„fc;,y)) 

y 

(18) 

where the initial and final state momenta k^ are cho- 
sen as independent variables, requires the interpolation 
in the momentum ky for the quantities on both sides of 
P34 that are calculated on the mesh {fci^}. The depen- 
dence on the discrete quantum numbers is suppressed 
since it is irrelevant for the consideration as well as the 
explicit form of function P34(fcz, fc^, y). k'y{kz,k'z,y) and 
and ky(kz, k'z, y) are the initial and final state Jacobi mo- 
menta ky expressed via fc^, fc^, and the angle between 
them y = • kz. We use the spline interpolation again 
with the spline functions Si^{k) (26l[27l[28i] such that for 
the function f{k), given on the mesh {ki^}, the values at 
any k may be obtained by 



/(fc)«^/(fc,J5.„(fc). 



(19) 



For P34 acting on the vector \Y) we obtain the following 
result (as a distribution) 



{k,kykz\P3i\Y) = J2 



S(ky -k,)r°°,2 , 

^ ^ 1 k'z^dk'z / dy 



i.2 



X Si^{ky{kz,k'z,y)) P34{kz,k'^,y) 



Ov 



E 



6{ky kj^ ) ^ 
kl 



(20) 



such that in the next step of the calculation, where the 
{kxkykz\P3i\Y) has to be multiplied by a smooth function 
f{ky) and integrated over ky, the result simply is the sum 



over the meshpoints {ki^} for the involved quantities, 

POO 

/ kldky f{ky){kxkykz\P3A\Y) ^y^f{k,^)Y,^{kx,kz). 

° iy 

(21) 

The integrations in Eq^jpO]) are performed using Gaus- 
sian integration rules [2^. The bound state pole (fTT]) 
is treated by the subtraction technigue much like the 
deuteron pole in the n-d scattering [241 • Note that the 
representation of the operators Pi and P34 is different 
from the one used in Refs. [l^, [2^ where final state mo- 
menta ky and kz were chosen as independent variables. 



III. RESULTS 

In order to calibrate our work we start by reproducing 
results of previous calculations, in particular the binding 
energy of ''He obtained with different realistic 2N inter- 
actions by different groups [13, [H, [13, [sH as well as the 
n-'^H phase shifts obtained with Mafliet-Tjon potential 
by the Grenoble group [s^l- Furthermore, we check the 
numerical stability of our calculations. These results are 
presented in the Appendix and show that the present 
algorithm is numerically highly reliable and capable of 
reproducing previous published results. 

Next we study the convergence of our calculations in 
terms of number of 2N , 3N, and AN partial waves us- 
ing the AV18 potential |33|] for the 2A^ interaction. In 
the calculations presented here for the n-^H scattering 
we include only the total isospin T — I states, but, 
within T = 1 , we take into account all couplings resulting 
from the charge dependence of the interaction. Including 
T — 2 states would yield an effect that is of 2nd order 
in the charge dependence and, therefore, is expected to 
be extremely small much like the effect of the total 3iV 
isospin T — ^ states in elastic n-d scattering. Coupling 
to T = 2 states is neglected also in all previous calcula- 
tions, but in configuration-space treatments the isospin 
averaging within T = 1 states is performed for the po- 
tential, whereas we perform it for the t-matrix. 

In Table H] we show n-^H phase shifts, 1~ mixing pa- 
rameter e, and total cross section at at En = 4 MeV 
neutron lab energy for increasing number of 2N partial 
waves. In all calculations we keep ly,lz < ^ and J < |. 
We apply additional restrictions that are different for 1+3 
and 2+2 states. We include all 1+3 states with Ix+ly < 8 
plus the states coupled to them by the tensor force; the 
above restriction is not applied if / < 2. We include 
all 2 + 2 states with Ix + ly + Iz < 10 plus states cou- 
pled to them by the tensor force. One finds that at least 
/ < 3 + '^P4 is needed for a well converged calculation. 
Likewise in Table|lT]we show similar results for increasing 
ly, Iz keeping / < 4 + ^G^ and J < |. At least ly, Iz < 
is needed to get quite satisfactorily converged results for 
the P-wave phase shifts, particularly ^P2. Finally in Ta- 
ble lllll we show results for increasing J, keeping ly, Iz < A 
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0~ (^Pn) 


1+ C^.Si ) 
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1" (^Pi) 


J- V -ri) 
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A A 00 
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I <2 + '^D-i 


-70.02 


24.43 


-62.05 


43.62 


22.65 


-UA6 


47.06 


2.301 


I <S + ^F4 


-69.68 


23.52 


-61.74 


43.37 


22.35 


-44.71 


46.71 


2.277 


7 < 4 + ^Gs 


-69.63 


23.62 


-61.69 


43.54 


22.38 


-44.69 


47.03 


2.288 


7 < 5 + ^ii-fi 


-69.61 


23.56 


-61.68 


43.53 


22.37 


-44.73 


47.00 


2.286 



TABLE I: n-"^H phase shifts, 1 mixing parameter e (in degrees) and the total cross section at (in barns) at En = 4 MeV 
neutron lab energy for increasing number of 2N partial waves and fixed ly,h < 4, J < |. The 2N potential is AV18. 





0+ ('5o) 


0- CPo) 


1+ es^) 


1- (^A) 


1- CP,) 


1- (^) 


2- (3P2) 


o-t 


ly,h < 


-69.70 




-63.50 










0.950 


^yi^z ^ 1 


-69.59 


22.62 


-61.94 


41.33 


22.65 


-44.49 


43.35 


2.163 


ly J Iz ^ 2 


-69.67 


23.19 


-61.75 


42.65 


22.05 


-44.88 


44.10 


2.196 


/ / < S 


-69.62 


23.65 


-61.72 


43.35 


22.34 


-44.84 


46.82 


2.279 


< 4 


-69.63 


23.62 


-61.69 


43.54 


22.38 


-44.69 


47.03 


2.288 



TABLE IL Same as in Table|T]for increasing ly,lz and fixed 7 < 4 -f "^Gs and J < |. The restriction on ly is not applied to 
3N partial waves with total angular momentum and parity J'" = 





0+ CSo) 


0- ePo) 


1+ Cs,) 


1- (^A) 


1- CP,) 


1~ W 


2- eP2) 


at 




-69.84 


23.95 


-53.98 


27.53 


17.55 


-9.48 


17.56 


1.268 


J<1 


-69.61 


23.26 


-62.41 


43.05 


22.34 


-44.85 


21.97 


1.715 


J<l 


-69.63 


23.61 


-61.69 


43.49 


22.37 


-44.63 


46.97 


2.285 




-69.63 


23.61 


-61.69 


43.53 


22.38 


-44.68 


46.99 


2.287 


- 2 


-69.63 


23.62 


-61.69 


43.54 


22.38 


-44.69 


47.03 


2.288 




TABLE 


III: Same as 


in Table H] for increasing J and 


fixed 7 < 4 + 


^Gs and ly, 


h < 4. 






0+ CSo) 


0- (^Po) 


1+ (^Si) 


1- (3Pi) 


1- CPi) 


I- (^) 


2- CP2) 


at 


AV18 


-66.12 


20.75 


-58.48 


40.09 


20.73 


-44.50 


42.51 


2.331 


Ref. [llj 


-66.5 


20.9 


-58.5 


37.3 


20.7 


-43.5 


41.0 


2.24 


Ref. [llj 


-66.3 


20.6 


-58.7 


38.6 


20.5 


-45.5 


40.1 


2.24 


rank 1 


-66.06 


26.99 


-58.55 


42.36 


22.15 


-44.81 


45.06 


2.488 


CD Bonn 


-64.63 


18.97 


-57.40 


39.44 


20.20 


-44.94 


42.47 


2.283 


Nijmegen I 


-65.61 


19.64 


-58.16 


39.62 


20.40 


-44.91 


42.13 


2.297 


Nijmegen II 


-65.98 


20.02 


-58.42 


39.69 


20.44 


-44.71 


42.22 


2.308 


INOY04 


-62.91 


16.73 


-56.00 


38.75 


19.47 


-44.55 


42.13 


2.216 


N3L0 


-65.54 


20.31 


-57.99 


40.94 


20.74 


-44.71 


43.98 


2.377 



TABLE IV: n-^U phase shifts, mixing parameter e, and total cross section at for AV18, CD-Bonn, Nijmegen I, Nijmegen II, 
INOY04, and N3LO potentials at En = 3.5 MeV together with resuhs from other calculations for AV18. We include 7 < 4-|-^G5, 
ly,l^ < 4, and J < |. 



and / < 4 -I- ^65. We find that the inclusion of at least 
J = I 3N states is necessary without which Ay has the 
wrong sign. Compared with previous calculations the 
present work exceeds in the number of 2N, 3iV, and 4iV 
partial waves included, providing very accurate results 



for all observables. 

In Table IIVI we show the results of the other calcu- 
lations for AV18 at En = 3.5 MeV which were com- 
piled in Ref. [iTJ. The present calculation confirms the 
work of the Grenoble and Pisa groups (second and third 
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lines, respectively) and clearly shows in the fourth line 
the shortcomings of the rank one representation of real- 
istic interactions calculated again using the present nu- 
merical algorithm. As in the work of Ref. 4] the total 
cross section gets to be at = 2A9 b which even slightly 
overestimates the experimental value. Calculations with 
other potentials, i.e., charge-dependent (CD) Bonn [sj], 
Nijmegen I and II [35|, inside- nonlocal outside- Yukawa 
(INOY04) potential by Doleschall [H ill, and the one 
derived from chiral perturbation theory at next-to-next- 
to-next-to-leading order (N3L0) [37|, show similar re- 
sults for all phases although N3L0 gives the largest P- 
wave phases leading to at = 2.38 b, the closest to the ex- 
perimental value at the resonance peak using two-body 
interactions alone. 

In Fig. [2] we show the total cross section for n-^K 
scattering as a function of the neutron lab energy; for 
clarity we skip the Nijmegen I and II predictions since 
they are between AV18 and CD Bonn. In the resonance 
region all potentials fail to reproduce the experimental 
data though some do better than others. As pointed out 
in Ref. 10] the nonlocal potential INOY04 that, by it- 
self alone, leads to the correct triton binding energy and 
slightly overbinds the a particle, shows the lowest total 
cross section at the peak. On the contrary CD Bonn 
and AVIS show higher total cross sections but also lower 
triton and a particle binding energies. 

In Table |V] we give the values for the triton and a 
particle binding energies, singlet and triplet scattering 
lengths oo and oi, and total cross section at En = 
and 3.5 MeV. The results we get for ao and ai agree 
with previous work for AV18 [13, El, and as shown in 
Fig.[3]correlate with the triton binding energy. Therefore 
interactions that lead to lower triton binding show the 
highest values for oq and ai and consequently the higher 
total cross sections at threshold. Nevertheless at En — 
3.5 MeV this correlation gets destroyed as the behavior of 
N3L0 shows. Further studies are needed to understand 
the features of N3L0 that give rise to this breaking of 
the correlation near the peak of the resonance. 



2.5 





£t 




ao 


ai 


at (0) 


at (3.5) 


AV18 


7.621 


24.24 


4.28 


3.71 


1.88 


2.33 


Nijmegen II 


7.653 


24.50 


4.27 


3.71 


1.87 


2.31 


Nijmegen I 


7.734 


24.94 


4.25 


3.69 


1.85 


2.30 


N3L0 


7.854 


25.38 


4.23 


3.67 


1.83 


2.38 


CD Bonn 


7.998 


26.11 


4.17 


3.63 


1.79 


2.28 


INOY04 


8.493 


29.11 


4.02 


3.51 


1.67 


2.22 



TABLE V: and ''He binding energies et and Sa (in MeV), 
n-^H scattering lengths ao and ai (in fm), and n-'^H total 
cross section at (in barns) at E„ — and 3.5 MeV neutron 
lab energy for different 2N potentials. 

In Fig. [3] we show the differential cross section da/dil 
and the neutron analyzing power Ay for n-'^H scattering 
at neutron lab energies of 1, 2, 3.5, and 6 MeV. In order to 



2.0 



1.5 




N3L0 
INOY04 
AVI 8 
CD Bonn 



0.1 



1.0 
En (MeV) 



FIG. 2: (Color online) Total cross section for n-^H scattering 
as function of neutron lab energy calculated with CD Bonn 
(solid curve), AV18 (dashed curve), INOY04 (dash-dotted 
curve), and N3LO (dotted curve) potentials. Experimental 
data are from Ref. (7|. 



5.0 



4.0 



J=0" 



AVI 8 N3L0 
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FIG. 3: Correlation between H binding energy et and n- H 
scattering lengths ao and ai. The predictions for AV18 (di- 
amonds), N3L0 (triangles), CD Bonn (circles), and INOY04 
(squares) are shown. The experimental data are from Ref. [38l ] 
(open circles), Ref. [3^ (open triangles), and Ref. (open 
squares) . 



get fully converged results we take into account all n-^K 
channel states with orbital angular momentum L < 3. 
The predictions of the four potentials differ mostly at 
forward and backward angles for the differential cross 
section and around the peak for the analyzing power. It 
is not obvious to us that the disagreement with the to- 
tal cross section data shown in Fig. [5] is compatible with 
the discrepancies we observe relative to the differential 
cross section data. Therefore it would be recommended 
that some of the experiments be repeated at specific en- 
ergies and Ay measured in order to further understand 
the implications of the 2N force models. 

One important observation that comes out of these cal- 
culations is the increased sensitivity of 4A^ observables to 
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changes in the 2N interaction. The variations due to the 
2N potential at the maximum of n-^Yi Ay lead to about 
16% fluctuations which are larger than the 10% fluctua- 
tions observed at the peak of Ay in low energy n-d scat- 
tering. This indicates that the AN system is more sensi- 
tive to off-shell differences of the 2N force than the 3A^ 
system. 

Finally, in Table IVII we investigate the possible cor- 
relations between the A^-puzzle in low energy n-d scat- 
tering and the underestimation of ctj in rt-'^H scattering 
in the resonance region. The experimental data for Ay 
in n-d scattering can be accounted for by a calculation 
with modified interactions in 2N ^Pj waves [s^, |4^. We 
use two models. The first one, AV18', is taken from 
Ref. [i^l; it corresponds to the AV18 potential that in 
^Pj waves is multiplied by strength factors 0.96, 0.98, 
and 1.06 for / = 0, 1, and 2, respectively. The sec- 
ond one, INOY04', is taken from Ref. [361] and differs 
from INOY04 by "^P/ wave parameters. Although both 
modified potentials provide quite satisfactory description 
of vector analyzing powers in low energy n-d scattering, 
they are incompatible with present day 2N data basis. 
E.g., the x'^/datum values with respect to the pp data, 
estimated using the Nijmegen error matrix [43j . i.e., by 
comparing to the Nijmegen phase shifts rather than to 
data directly, are 3.5 for INOY04' and 4.4 for AV18' po- 
tentials. However, those modifications of the potentials 
are unable to resolve the at discrepancy in n-^H scatter- 
ing. The at is slightly increased for AV18' but it gets 
even lower for INOY04', indicating that at depends on 
the 2N ^Pj wave interaction in a different way than the 
Ay in the n-d scattering. 

IV. CONCLUSIONS 

In the present paper we developed a new numerical 
approach to solve four-nucleon scattering equations in 
momentum-space. The method uses no uncontrolled ap- 
proximations, is numerically very efficient and therefore 
can include very large number of partial waves, thereby 
yielding well converged and very precise results. The 
developed approach is applied to n-^B. scattering below 
three-body breakup threshold. The calculations with 
various realistic 2N potentials underestimate the total 
n-'^H cross section data in the resonance region as al- 
ready found by other groups. However, probably due to 
the inclusion of more partial waves, the numbers we get 
are slightly higher and closer to the data; they also de- 
pend on the choice of the 2N potential. The new results 
also show that 4N observables are more sensitive than 



3iV observables to the off-shell nature of the 2N inter- 
action. Furthermore, the modifications that are required 
to introduce at the level of the '^Pj 2N partial waves to 
remove the discrepancies in n-d Ay at low energy, do not 
remove the disagreement observed in the total n-^H cross 
section around En = 3.5 MeV. Finally, to understand the 
compatibility between existing n-'^H total and differential 
cross section data it would be advisable to repeat some 
of those experiments at specific energies. 
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APPENDIX A 

As mentioned in Section IIIII we present here our re- 
sults for the binding energy of ^He and n-^H phase shifts 
obtained with Mafiiet-Tjon potential as well as the nu- 
merical stability check of our results. In Table IVIII we 
show the a particle binding energy for increasing number 
of 2A^ partial waves and compare with previous works. 
Results with AV8' are calculated without the Coulomb 
interaction in order to compare with Ref. [3i3 |. On the 
other hand calculations from Ref. [l^, |3l| with CD Bonn 
include coupling between total isospin T = 0, 1 and 2 
states while we consider only T = 0. In contrast to our 
scattering calculations, here we perform the isospin av- 
eraging not for the t-matrix but for t he p otential like it 
has been done in calculations of Ref. [10l| . Overall these 
results indicate that our algorithm is accurate and reli- 
able. 

For n-^H scattering we compare in Table IVIIII the 
results of our calculations with those of the Grenoble 
group f33|. Again our phase shifts agree within a few 
tenth of a degree or better leading to identical total cross 
sections over a wide range of energies. 

In Table HXl we demonstrate the stability of our results 
increasing the number Nk of momentum meshpoints. All 
calculations use AVIS potential with 27V partial waves 
/ < 4 ^Gs, Zy, ^2 < 4 and J < |. The results with 
A^fc — 20 which is the choice for calculations of Tables U 
— nil] are converged to better than 0.05 deg for all phase 
shifts. 
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INOY04 
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INOY04' 
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19.06 
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42.21 
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TABLE VI: n-^H phase shifts, 1" mixing parameter e, and total cross section at at En = 3.5 MeV for original AV18 and 
INOY04 potentials and their versions AV18' and INOY04' with modified P-waves. 
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TABLE VII: ■'He Binding Energy (MeV) for increasing number of 27V partial waves 


characterized by maximal total angular 


momentum 7. 
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TABLE VIIL n-^H phase 


shifts at different neutron lab energies for the Mafliet-Tjon potential. 




Nk 0+ CSo) 


0- CPo) 


1+ (^Si) 


1 


- 


) 


1- (iPi) 


1- (e) 2- CP,) 


o-t 


14 -70.03 


23.78 


-62.01 




43.49 




22.35 


-44.58 46.93 


2.290 


15 -69.57 


23.57 


-61.61 




43.64 




22.40 


-44.72 47.18 


2.292 


16 -69.43 


23.51 


-61.51 




43.63 




22.38 


-44.75 47.19 


2.290 


18 -69.72 


23.66 


-61.76 




43.55 




22.39 


-44.67 47.03 


2.289 


20 -69.63 


23.62 


-61.69 




43.54 




22.38 


-44.69 47.03 


2.288 


22 -69.68 


23.64 


-61.73 




43.53 




22.39 


-44.68 47.01 


2.288 


24 -69.66 


23.64 


-61.73 




43.53 




22.39 


-44.68 47.00 


2.288 


25 -69.67 


23.64 


-61.73 




43.53 




22.39 


-44.68 47.00 


2.288 



TABLE IX: n-^H phase shifts, 1" mixing parameter e, and the total cross section at at En = 4 MeV for increasing number of 



meshpoints Nk for the momenta 



The 2N potential is AV18. 
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